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Spontaneous breaking of translational symmetry—known as ‘density wave’ order—is common in 
nature. However such states are strongly sensitive to impurities or other forms of frozen disorder 
leading to fascinating glassy phenomena. We analyze impurity effects on a particularly ubiquitous 
form of broken translation symmetry in solids: a Spin Density Wave (SDW) with spatially modulated 
magnetic order. Related phenomena occur in Pair Density Wave (PDW) superconductors where 
the superconducting order is spatially modulated. For weak disorder, we find that the SDW / 
PDW order can generically give way to a SDW / PDW glass—new phases of matter with a number 
of striking properties, which we introduce and characterize here. In particular, they exhibit an 
interesting combination of conventional (symmetry-breaking) and spin glass (Edwards-Anderson) 
order. This is reflected in the dynamic response of such a system, which—as expected for a glass— 
is extremely slow in certain variables, but—surprisingly—is fast in others. Our results apply to all 
uniaxial metallic SDW systems where the ordering vector is incommensurate with the crystalline 
lattice. In addition, the possibility of a PDW glass has important consequences for some recent 
theoretical and experimental work on La 2 ~ x Ba x Cu 2 04 . 


I. INTRODUCTION 

A variety of electronic solids settle into equilibrium 
states that spontaneously break the translational sym¬ 
metry of the underlying crystal*. Well known examples 
are Charge Density Wave (CDW) and Spin Density Wave 
(SDW) orders where either the electron’s charge or its 
spin forms a frozen periodically oscillating pattern. Den¬ 
sity wave order has been found in conventional metals 
as well as in strongly correlated systems such as the un¬ 
derdoped cuprates, iron pnictides, and organic materials, 
and are intertwined with many other fascinating phenom¬ 
ena such as for instance high temperature superconduc¬ 
tivity. 

It has long been recognized that density wave orders 
of various kinds are strongly sensitive to the presence of 
impurities. There is a large literature on the fascinat¬ 
ing effects of quenched disorder on charge density wave 
systems. In contrast, despite the common occurrence of 
spin density wave ordering, surprisingly little attention 
has been devoted to impurity effects on SDW systems, 
and this is the subject of this paper. 

It is important to distinguish between collinear SDWs 
where the spin orientation oscillates in space along a fixed 
common direction and spiral SDWs where the spin ro¬ 
tates around an axis as a function of space while the 
magnitude |5(r)| is constant. Both kinds of SDW order 
break both spin-rotation and lattice translation symme¬ 
tries, but the latter retains a combination of the two as a 
symmetry. As a consequence, disorder effects on collinear 
SDWs are stronger and are the focus of our study. We 
show that weak non-magnetic disorder transforms the 
SDW state into a new glassy state of matter—distinct 
from the conventional spin glass—which we dub the SDW 
glass (see Fig. [I]). In addition to the characteristics of a 


conventional spin glass, i.e. the presence of locally frozen 
moments but absence of long-range spin order, the SDW 
glass spontaneously breaks spin-rotation symmetry and 
hosts an associated Goldstone mode (see Fig. [2]). 



FIG. 1: Magnetic moments in the ground state of the SDW 
glass. The disorder pins the domain walls into a random con¬ 
figuration, but the structure of ‘anti-phase’ domain walls per¬ 
sists. Upon crossing of each domain wall, the local magne¬ 
tization changes sign. Thus the disordered state inherits the 
collinear structure from the parent SDW, where the axis along 
which the moments point is selected sponteneously, breaking 
the spin rotation symmetry. 

Questions closely related to the ones we study arise 
in considering the effects of impurities on Pair Den¬ 
sity Wave (PDW) superconductors of the kind proposed 
to be realized in the high temperature superconductor 
La 2 - x Ba x Cu 2 0 A (see Ref. [2]for a recent review). These 
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FIG. 2: Goldstone mode in the SDW glass. The configuration 
of the domain walls is identical to the ground state, but the 
orientation of the moments slowly varies in space 


are superconducting states where the pair amplitude is 
modulated in space. Such states were originally con¬ 
ceived by Larkin and Ovchinnikov 3 when considering the 
effects of a Zeeman magnetic field on an s-wave supercon¬ 
ductor (closely related states were considered by Fulde 
and FerrelP) . Motivated by the phenomenology observed 
in La 2 - x Ba x Cu, 204 , Refs. 06] proposed that a PDW 
state is realized at zero magnetic field in this system. Its 
microscopic origin remains to be understood. Like the 
analogous SDW order, the PDW state is also expected 
to be strongly sensitive to impurities. Refs. l5|6l sug¬ 
gested that disorder necessarily introduces half-quantum 
(he/ 4e) vortices leading to broken time reversal symme¬ 
try. We revisit this issue and find that in the phase anal¬ 
ogous to the SDW glass, such vortices will not be induced 
by weak disorder and that time reversal symmetry is pre¬ 
served. A fluctuating version of the PDW state also ap¬ 
pears as a ‘mother’ state that controls the physics of the 
pseudogap regime in the cuprates in a recently proposed 
theory 7 . 


II. SDW AND PDW: ORDER PARAMETERS 
AND TOPOLOGICAL DEFECTS 

For simplicity we focus initially on a unidirectional 
collinear SDW at wave-vector Q in which the spin os¬ 
cillates as 


S(r) = So cos (Q ■ r) N, (1) 

where A is a real unit vector and So is the amplitude. 
We further assume that Q is incommensurate with the 
crystalline lattice. Such SDW order is sometimes also 
referred to as a ‘spin stripe’. A simple Landau argument 
shows that this pattern of spin ordering induces CDW 


order at wave vector 2 Q: 

Pgdw ~ P2Q cos (2 Q ■ r ) (2) 

To discuss situations in which the SDW order is fluc¬ 
tuating (either due to quenched disorder or due to ther¬ 
mal/quantum fluctuations) we write the spin as 

S(r) ss Sb cos(Q • r + 9(r))N. (3) 

9 describes the phase of the SDW (or the stripe displace¬ 
ment, in the stripe picture). We take both 9 and N 
to be space (and possibly time) dependent but varying 
on length scales long compared to the SDW wavelength. 
The fluctuating SDW order parameter near wave vector 
±Q is thus 


S±q(x) = S 0 e ± ^N(x). (4) 


The corresponding fluctuating CDW order parameter is 


_ _ _0 c ±2i0 

P±2Q ~ P2Q e 


( 5 ) 


It is useful to formulate discussions of fluctuations in 
terms of the two separate fields b = e™ and N. Both 
the SDW order parameter Sq ~ bN and the CDW or¬ 
der parameter p 2 Q ~ b 2 are composites made out of b 
and N. Clearly the b , N representation has a Z 2 gauge 
redundancy under b -A — b.N -A —IV so that neither of 
them are directly physicaP®!. The SDW/CDW order 
parameters are of course gauge invariant. 

Let us now turn to the closely analogous PDW 
state (sometimes called a striped superconductor or the 
Larkin-Ovchinnikov phase). This is a superconductor in 
which the pair wave function A is modulated in space: 


A(r) = A 0 cos(Q ■ r) (6) 

This too will induce CDW order at wave vector 2 Q. 
When fluctuating we may write 

A(r) « A 0 cos(Q • r + 0(r))e l< ^ r) (7) 

The Fourier components of A(r) near ±Q are thus ~ 
e *O±0). There is again a Z 2 gauge redundancy under 
9 —> 9 + ir, <j> —> cj) + n. 

The PDW state is thus conceptually very similar to a 
SDW state with just XY spin symmetry. However there 
is an important difference in the action of time reversal 
symmetry. The PDW state preserves time reversal while 
the SDW breaks it. Formally this is because the C/(l) 
charge conservation symmetry (broken in the SC) does 
not commute with time reversal while spin rotations do. 
Nevertheless we will consider both orders within the same 
framework. Unless otherwise specified we will phrase our 
discussion in terms of SDW order. 

The structure of topological defects in these density 
wave stateP^ 2 will play a crucial role below. Of partic¬ 
ular importance are dislocations in the CDW pattern. 
These are line defects in 3D and point defects in 2D 
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around which 9 winds. The ‘elementary’ strength-1 dis¬ 
location where the CDW phase 29 winds by 2n requires 
that N winds by 7r so that the SDW order e l0 N is single 
valued. In contrast strength-2 dislocations have 29 wind 
by 47 t without any winding of N. 

Exactly the same considerations also apply in the su¬ 
perconducting context as described in Refs. 131131141 It is 
interesting to consider the physical interpretation of the 
various topological defects in this case. The strength-1 
CDW dislocation now requires that the superconducting 
phase 4> wind by n. This corresponds to a superconduct¬ 
ing vortex with magnetic flux i.e half the usual flux 
quantum. Strength-2 CDW dislocations in contrast do 
not bind to superconducting vortices. 


III. IMPURITIES: MODELS AND 
PRELIMINARIES 

We want to consider the fate of the SDW in the pres¬ 
ence of weak non-magnetic impurities. Such impurities 
lead to a random potential V ( x ) that couples linearly to 
the CDW order parameter, i.e. 

Hdis. = F(x) ■ V9 + V 2 * Q (x)e 2ie + V 2Q (x)e~ 2ie , (8) 

The first term couples to the long wavelength part of 
the charge density (with F random) and the second 
to the density near the ordering wave vector. Here 
V 2 q(x ) = /q~ 2 Q eiq ^(q)- There is however no linear 
coupling to the primary SDW order parameter. The im¬ 
purity coupling is captured by a simple lattice model: 

H = - \ E % • S Qj e^ - V ~Y. + c ' c 

<ij> i 

= - J Ni ■ Nj cos (Oi - 9j + rjij) 

<ij> 

— u^cos(2 9i — ai) (9) 

i 

Here oti,rjij are random uncorrelated variables. Other 
equivalent lattice models may be formulated and are de¬ 
scribed in the Appendix [B] 

As is well known 15 16 the ‘random field’ disorder de¬ 
stroys LRO in the CDW for physical dimensions d < 4. 
The elastic energy cost of adjusting to disorder over a 
scale L scales as L d ~ 2 while the energy gain due to the 
disorder potential scales as L d / 2 , thus the latter dom¬ 
inates for d < 4. As a consequence, beyond a length 
scale (known as the Larkin length) ~ {J/v) 2 ^ A ~ d ' > 
long range CDW order is destroyed. This immediately 
implies the absence of long range SDW order as well 
(as long range SDW order if present would have induced 
CDW order). At distances longer than the disordered 
SDW enters a phase of matter that we dub the SDW 
glass and whose physics we describe below. 


IV. SDW GLASSES IN 3D 

We begin our analysis in d = 3 dimensions by reviewing 
the physics of pinned CDW systems. In pioneering work, 
Ref. [TTI proposed that at weak disorder the pinned CDW 
enters an ‘elastic glass’ phase where long dislocation loops 
do not occur. This has been substantiated by numer¬ 
ical calculation:^ anc j by general scaling argument^^. 
The resulting state is described by a random field XY 
model for the CDW order parameter where dislocations 
are suppressed. Many approximate treatments, notably a 
sophisticated Functional Renormalization Group (FRG) 
calculation 1 ®^, show that the CDW order parameter de¬ 
velops power law correlations: 

P*2q{x) P2 q{x') ~ ^ ^ ( 10 ) 

The exponent d c is universal. To leading order in the e 
2 

expansion, d c = ^. Thus for d = 3, d c is estimated 
to be ~ 1.1. The power law decay of the spatial CDW 
correlations implies power-law Bragg peaks in the CDW 
structure factor 

'S'cdw( 2<5 + <5<j) ~ \5q\ da ~ 3 , (11) 

The elastic glass phase is therefore also known as the 
‘Bragg glass’ phase. 

Let us now consider the implications for the SDW or¬ 
der. The impurities do not couple linearly to N but will 
lead to random exchange energies (‘random bond disor¬ 
der’). However the absence of dislocations in 9 means 
that there is no frustration of the collinear ordering of 
N. Thus at weak disorder N will continue to have true 
long range order. The SDW order will has power law 
correlations inherited from the correlations of e lS . 

5q(x) • Sq(x') = e* 6l ( x )e _i6l ( x ')JV(x) • 7V(x') 

^ gi0(x)g — 20(x') 

1 

|x — x'l^ 

The exponent d s can be estimated within the FRG in d = 
4 — e. At leading order in the e expansion we have d s = 

2 i 

^. This is j of the CDW exponent d c , as the probability 
distribution for 9 is Gaussian to this order (for a recent 
discussion see Ref. I2B- This gives the estimate d s ~ 0.27 
in d = 3. However beyond leading order the distribution 
will not be Gaussiarpll, and hence in general d s ^ c ^f. 
Note that the SDW correlations decay much slower than 
the CDW correlations. This is expected since the CDW 
order is the one directly affected by the disorder. 

Correspondingly, the spin structure factor exhibits 
power law Bragg peaks 

S SDW (Q + 5q)~\5q\ d °- 3 . (12) 

This power law Bragg peak should be visible in neutron 
diffraction measurements on weakly disordered SDW ma¬ 
terials. 
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But what does it mean for N to be ordered? N 
is not gauge invariant and hence not directly observ¬ 
able. However ordering of N implies ordering of the spin 
quadrupole moment 


Q af 3 = N a N 0 - 6 -^N 2 (13) 

Thus even though the SDW order is destroyed long range 
spin quadrupole order (also known as spin nematic or¬ 
der) is preserved. The system develops spontaneous spin 
anisotropy without long range SDW ordering. 

This SDW glass phase has a simple physical descrip¬ 
tion. The spins are frozen in time but the phase of the 
SDW is randomly disordered in space. The spin nematic 
order means that the spins retain a common axis along 
which they randomly point up or down. The freezing 
of the spins means that there is a non-zero Edwards- 
Anderson spin glass order parameter 


9ea = lim (S(x, 0) • S(x, t)) (14) 

£—>■00 

= lim (cos (0(x, t) — 6{x, 0)))(iV) 2 ^ 0, (15) 

£—>•00 

The disordered SDW is thus a uni-axial spin glass in a 
Heisenberg spin system with the axis of spin orientation 
determined spontaneously. It is clearly distinct from the 
conventional Heisenberg spin glass. 

The spin nematic order in the SDW glass phase leads to 
propagating Goldstone modes (nematic director waves). 
The structure of magnetic moments characteristic for 
the ground state and for a soft excitation are shown in 
Figs. [I |2| This should be contrasted with the Halperin- 
Saslow22T spin wave modes in a Heisenberg spin glass 
which are typically damped. 



FIG. 3: In the SDW Bragg glass, the anti-phase domain walls 
are pinned in a configuration without dislocations. 

The SDW glass exhibits striking differences from con¬ 
ventional spin glasses when placed in a weak magnetic 
field. Specifically, we will contrast the SDW glass to a 


Heisenberg spin glass (HSG) given by 

#HSG = Jx,yS x ■ Sy , (16) 

xy 

where J x y are random exchange couplings. Recall that 
in the absence of disorder, a collinear SDW aligns itself 
perpendicular to the applied field, with a slight canting 
of the moments, like an antiferromagnet. In the isotropic 
SDW glass, the canting is achieved by a rotation of 
N. This corresponds to a fast (Goldstone) mode and 
therefore does not exhibit slow (glassy) dynamics. In 
the HSG such a fast mode does not exist and the time 
scales for adjusting to the Held are necessarily long. 

For an anisotropic SDW glass N remains pinned to a 
specific direction in a weak enough magnetic field, and 
the system can only respond via the slow dynamics 
of e . The cross-over value of the magnetic field B c 
between fast and slow dynamics is thus given by the 
strength of the anisotropy. In the HSG, a similar 
cross-over occurs, but at a much larger scale, determined 
by the typical exchange coupling J. 

The same analysis presented here for the SDW also 
applies for a PDW. The analog of the spin nematic order 
parameter is a uniform charge-4e superconducting order 
A 4e . However we are not aware of any system that is 
proposed to host a PDW at zero magnetic field in d = 3 
dimensions. Our results should be pertinent though to 
Larkin-Ovchinikov pairing induced by a magnetic field in 
three dimensional superconductors. The breaking of time 
reversal allows additional terms in the Hamiltonian, in 
particular a linear coupling between the gradients of the 
phases of CDW and the condensate, i.e. in the continuum 
limit 

H b = A f V0 • V6» (17) 

J X 

However, the PDW glass phase is perturbatively stable 
against such a term. Clearly the ground state in the 
absence of disorder is unaffected by Hb for small A. In 
the presence of disorder, at long distances L £ L , H B 
contributes to the random bond energy for the charge-4e 
superconducting order parameter A 4e , which is irrelevant 
in the PDW glass where (A 4e ) ^ 0. The same conclusion 
can be readily obtained through an FRG analysis (see 
Appendix [A} . 


V. SDW AND PDW GLASSES IN 2D 

We now turn to d = 2. Once again the random field 
will destroy long range CDW and hence long range SDW 
order. The fate of dislocations is however more subtle. 
In the simpler problem of the 2D random field XY model 
(appropriate to describe unidirectional CDW ordering 
not derived from a more primary SDW or PDW order), 
topological defects always proliferate at long scales, lead¬ 
ing to exponentially decaying correlations. These defects 
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cost elastic energy which must be balanced against the 
energy gain due to the random correlated potential in¬ 
duced by the random field. Refs. !23l2il show that at long 
enough scales the optimized potential energy for intro¬ 
ducing vortices dominates so that it is always favorable 
to nucleate defects. 

In the SDW system, single and doubled CDW disloca¬ 
tions have different elastic cost—the energy of the former 
depends on the stiffness associated with spin distortions 
while the energy of the latter does not. Hence they could 
potentially behave very differently. While one expects 
that doubled dislocations are always generated at long 
length scales the fate of isolated single dislocations is less 
clear. 

Let us first describe a putative state where doubled 
dislocations have proliferated but single ones have not. In 
such a state the CDW correlations, and hence the SDW 
correlations decay exponentially. Despite that there is 
long range spin nematic order. Thus this is a 2D SDW 
glass phase with coexisting spin nematic order. 

It is particularly interesting to note the meaning of 
these issues in the PDW context. The analog of the 
spin nematic then is a uniform charge-4e superconductor. 
Further (see Section 0 a strength-1 CDW dislocation is 
bound to a 7T winding of the SC phase. This corresponds 
to a half-quantum vortex with flux. Clearly two dif¬ 
ferent dislocations of this kind are possible depending on 
whether the flux is positive or negative. Any such sin¬ 
gle dislocation necessarily breaks time reversal symmetry. 
Thus if disorder nucleates these single dislocations then 
the PDW glass will break time reversal spontaneously by 
generating randomly placed ±^| half-quantum vortices 
of either sign. Such vortices can be imaged using local 
probes of magnetism such as a scanning SQUID micro¬ 
scope, and can serve as a key experimental test of the 
proposed PDW state in La 2 - x Ba x Cu 204 : . 

We therefore now pose the question of whether such 
single dislocations are necessarily generated at weak dis¬ 
order in either the SDW or PDW system. 

For simplicity and to provide a unified discussion of 
both SDW and PDW systems, we will specialize to XY 
spins. The fate of dislocations may be discussed within 
an appropriate elastic model which takes the form 

H=J d 2 x^(V(t>) 2 + ^{V6- f) 2 -vcos{26~a) 

(18) 


Here K c , K s are the stiffnesses of the CDW and SDW 
(or PDW) order parameters respectively. The phase of 
the SDW (or PDW) order parameter is ± 0, and that 
of the CDW order parameter is 26. The disorder is 
taken to be delta-correlated fi(x)fj(y) = DpSijS(x — y), 
a(x)a(y) = S(x — y). Eq. (18) is the basis for our anal¬ 
ysis of dislocations in the rest of this section. 

The 6 sector is again described by an XY model with 
random anisotropy. However, in this case there is no sta¬ 
ble dislocation free phase even for weak disorder. We 
will begin with discussions by reviewing some basic facts 


on this model: In the case of random forces (i.e the / 
term), but without random fields, the relative displace¬ 
ment 9 X — 6 X ’ grows logarithmically with distance, with 
a coefficient given by the variance Dp of the random 
forcd^. The energy cost of the cheapest dislocation is— 
E ~ K c { 1 — y^Dp/Dp) log L, i.e. there is a critical 
strength Dp of the random force, below which no iso¬ 
lated dislocations are present in the ground state. In the 
presence of random fields but with dislocations excluded, 
Dp is renormalized without bouiicP^ as Dp = ClogL. 
The coefficient C is temperature dependent. This implies 

(9 X — 9x') 2 ~ log 2 \x — x'\. (19) 


Thus C n (r) = ^e m U x+r )e -m U x )) decays faster than 
any power law, and the associated Bragg peaks are 
indistinguishable from the case of short-range correla¬ 
tions. When dislocations are allowed, they always be¬ 
come relevanlP^il f or l ar g e enough L since Dp(L) ~ 
logL > Dp. Beyond the scale where dislocations pro¬ 
liferate, C n {r) decays exponentially. 

"kjk. 

Fully disordered 

- (all dislocations) 


V 

—> 


SDW or X 
PDW Glass \ 

(no single dislocations) 


FIG. 4: Schematic phase diagram for a disordered SDW (or 
PDW) in 2D as a function of disorder strength v and the ratio 
K a /K c of the stiffnesses associated with spin waves / phonons 
in the absence of disorder 

The analysis in Appendix [C| shows that single disloca¬ 
tions proliferate at a scale 

£p) ~ ^gVirod+^V^cPioga/a (20) 

Doubled dislocations on the other hand proliferate at a 
different length scale 

$) » lo «U/a ( 21 ) 

It follows that for large -jp strength-2 dislocations prolif¬ 
erate at a shorter length scale. In that case the physics 
for distances longer than Y v is modified from that de¬ 
scribed by the elastic model. The prevalence of doubled 
dislocations renormalizes K c to zero without affecting 
K s . The net energy cost of the spin distortion associ¬ 
ated with a strength-1 dislocation is then jK s ln(^) but 
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FIG. 5: Summary of numerical data. Left: Spatial decay of the correlations functions Cinematic and Ccdw- The inset shows 
a log — log plot of the same data. Combined error due to sampling and disorder averaging is below the size of the shown 
data points. The qualitative feature that nematic correlations decay much slower than CDW correlations is generic for low 
temperatures and weak disorder. Center: Temperature dependence of the spin-nematic correlation length ^Nematic shown for 
parameters corresponding to a SDW glass phase (green) or a fully disorderd phase (purple) at low temperatures. For the first 
set of parameters, ^Nematic grows rapidly as temperature is lowered, indicating a phase transition into the SDW glass phase. 
For the second set of parameters, £jv shows no sign of divergence and is expected to saturate at a finite value such that the low 
temperature glassy phase is smoothly connected to the high temperature phase. Solid lines are drawn as a guide to the eye. 
Right: The stiffness pN for spin nematic fluctuations, as a function of system size shown for both the SDW glass and the fully 
disordered phase. In the SDW glass (k c = —0.1k,,), pjv depends only weakly on system size L and always takes a finite value, 
while for k c = 0.1k s it exhibits substantial sample dependence and rapidly decays with system size. 


now the optimal potential energy gain of the dislocation 
will just be a constant (see Appendix 0- Thus in this 
case strength-1 dislocations are suppressed even at the 
longest length scales. In the SDW glass this implies that 
spin nematic long range order survives at weak disorder 
in this regime. In the PDW glass in this regime time 
reversal is preserved and the disorder does not introduce 
half-quantum vortices. Rather there is true long range 
order in the charge-4e superconducting order parameter. 

In the opposite regime of large Ip single dislocations 
will proliferate first even at weak disorder. In either 
regime at strong disorder single dislocations will prolifer¬ 
ate. The schematic phase diagram is shown in Fig. [4] 


A. Numerics 

To support the arguments presented above we per¬ 
formed Monte-Carlo simulations of the Hamiltonian 

HyiC — «c ^ ' COs(20i 2$j -(- ) 

(O') 

- K s ^2 COS (8i - 9j + 7)ij) COS (<j>i - (f)j) 

(o> 

+ u^cos (20i + Qi)). (22) 


where (r, r') denotes nearest neighbors on a square lat¬ 
tice, a{r) is a uniform random variable for each site and 
is a gaussian random variable for each link with vari¬ 
ance Dp. We note that rj r y is generated under renor¬ 
malization and does not need to be explicitly added to 
the Hamiltonian. However, it is convenient for study¬ 


ing the case of weak random field. There, the large dif¬ 
ference between and £y makes numerical simulations 
challenging. Taking Dp < Dp significantly reduces 
while otherwise retaining the same physics as the Dp = 0 
model. 

Configurations were generated using the standard Hy¬ 
brid Monte Carlo algorithrrPT The correlations functions 


C'cDw(bj) = (exp (2 i9i - 2 iOj)) (23) 

C]Sfematic(b j) = (exp 2^(/>j)) (24) 


were measured, as well as the nematic stiffness parameter 
Pn, i.e. the response of the free energy F to a twist $ in 
the boundary conditions for <f>: 


Pn 


1 


d 2 


k s + 4k c (94> 2 


Fm 


•5=0’ 


(25) 


where in the absence of disorder, lim p s = 1. pm was 

T-K) 

obtained directly as a non-local correlations function on 
systems with periodic boundary conditions. The data is 
summarized in Fig. [5] 

As expected, equilibration becomes difficult to achieve 
at low temperatures due to slow glassly dynamics, 
putting severe constraints on the achievable system sizes. 
We ensured that equilibration is indeed achieved by us¬ 
ing three different initializations and confirming that the 
measurements are independent of this choice. We took 

(i) disordered configuration with <p, 8 completely random, 

(ii) ordered configuration with (e l6 ) = (e®^) ^ 0, (iii) 
an ‘annealing’ protocol where the temperature is succes¬ 
sively lowered towards the target value. 
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a. Correlation functions - Typical results for the 
two correlations functions Ccdw and CNematic are shown 
in Fig. [5] (left) for (3n s = lOfiv = 5, n c = v measured for 
15 disorder realizations on a system of size 160 x 160. For 
these system sizes, the spatial correlations exhibit very 
little sample dependence, and the combined error due to 
sampling and disorder averaging is smaller than the sym¬ 
bols used to plot the data. More generally, Ccdw always 
exhibits rapid exponential decay at low temperatures and 
moderately weak disorder strength, while Carnatic decays 
much slower. Depending on the parameters Carnatic can 
also be exponential, or follow a power-law up to achiev¬ 
able system sizes. 

b. Temperature dependent correlation length At 
moderate temperatures, CNematic is always follows an 
exponential decay with a temperature dependent cor¬ 
relations length £at, shown in Fig. [ 5 ] (center) for both 
k c = — 0.1k; s (SDW glass) and for k c = 0.1k s (fully disor¬ 
dered). At high temperatures £/v(re c = 0.1 k s ) > £/v(k c = 
—0.1/c s ), as would be the case for v = 0 (no disorder). As 
temperature is lowered, disorder becomes important and 
affects the relative size of the Cv For parameters corre¬ 
sponding to a SDW glass phase at lower temperatures, £jv 
diverges. For parameters where single dislocations even¬ 
tually proliferate, the correlation length shows no sign of 
divergence. 

c. Helicity modulus The stiffness parameter pn as 
defined in Eq. (251 provides a sharp distinction between 
the two low temperature phase. A non-zero value for 
Pn in the thermodynamic limit is a characteristic fea¬ 
ture of the SDW glass phase while in the fully disordered 
phase where all vortices have proliferated lim pn = 0. 

L —>00 


Results for the stiffness are shown in Fig. [ 5 ] (right) for 
/3k s = 10/3v = 5, re c = ±v. In the SDW phase, pn ex¬ 
hibits little variation between different samples and ap¬ 
proaches a finite value at the largest system sized that we 
considered. In the phase where single dislocations pro¬ 
liferate, pn fluctuates strongly from sample to sample, 
around an average that rapidly decays with increasing 
system size. 

Despite the limitations on sample size, the result of 
our simulations are fully consistent with the analytical 
predictions. In particular we find that depending on 
our choice of parameters, the system exhibits qualita¬ 
tively different behavior, corresponding to two different 
phases at low temperatures—the SDW glass phase which 
we introduce here, and a more conventional glassy state 
where all correlations decay exponentially. We expect 
that these conclusions remain valid for L —> 00 and 
T-» 0. 


VI. DISCUSSION 

Our results should be pertinent to a wide variety of 
systems. Below we highlight a few specific interesting 
examples. 

SDW ordering is very common in electronic solids, and 


is often incommensurate with the underlying lattice. The 
classic example 28 is elemental Chromium Cr. There have 
been several studies of the suppression of SDW order in 
Cr when it is alloyed with other transition metal ele¬ 
ments (see Refs. 1291301 and references therein), for in¬ 
stance Vanadium V . At low V concentrations where the 
disorder is weak, we expect our results to apply directly 
and predict the occurrence of the 3D SDW glass phase. 
Consequently the static spin structure factor is a power 
law (Eqn. 12), potentially visible in high resolution neu¬ 
tron diffraction studies. 


The SDW glass offers an interesting experimental op¬ 
portunity to probe the physics of the 3D Bragg glass in 
magnetic systems. The original theoretical proposa l 17 * 31 ! 
of the Bragg glass phase spurred a search for it in a few 
experimental systems, notably in vortex matter inside su¬ 
perconductors (for a review see Ref. [32]). Ref. f33l provided 
evidence for the predicted power law Bragg peaks in a dis¬ 
ordered vortex lattice through small angle neutron scat¬ 
tering . Experimental evidence for Bragg glass physics in 
CDW systems seems scarce—probably due to the strong 
coupling to disorder of the CDW order parameter. Re¬ 
cently however Scanning Tunneling Microscopy images of 
the quasi-two dimensional CDW ordered system NbSe 2 
have been interpreted^ in terms of a Bragg glass picture 
which might describe intermediate length scale physics. 
SDW systems of the kind considered in this paper offer 
a different context for Bragg glass physics which may be 
more directly amenable to experimental studies. 

We emphasize that the SDW glass is distinct from the 
conventional Heisenberg spin glass. The SDW glass is 
also distinct from the ‘cluster spin glass’ which rnacro- 
scopically is the same phase as the usual spin glass. 
It is i ntere sting that even some ‘classic’ metallic spin 
glassed 35 3fi - (for instance CuMn) actually have substan¬ 
tial short range SDW orded^ (visible in neutron diffrac¬ 
tion as a well defined finite wave vector peak). Physically 
these are usefully understood as obtained from local pin¬ 
ning of SDW fluctuations of elemental Cu around Mn 
impurities. A useful theoretical approach to understand¬ 
ing these systems may be to start with the SDW glass 
described in this paper and then to disrupt it with topo¬ 
logical defects at long length scales. 

Turning to 2d systems, it is interesting to consider 
the very lightly doped cuprates within the framework 
of our results. At low-T, these have long been reported 
to have spin glass order but also show substantial co¬ 
existing SDW correlations 3 ^. As we have argued if the 
‘parent’ SDW order is uniaxial ( i.e not a spiral) then 
two distinct kinds of glassy states are possible in both of 
which the SDW correlations decay exponentially on long 
scales. When only doubled dislocations of the accompa¬ 
nying CDW are induced (the 2D SDW glass phase), spin 
nematic order persists. If however single dislocations are 
also induced the resulting phase is smoothly connected 
to the conventional spin glass. 

What experiments can help distinguish between these 
two phases? As we discussed they will have rather sim- 










ilar peaks in neutron diffraction. More telling will be 
local probes of the spin dynamics for instance through 
NMR. The local dynamic spin susceptibility in the SDW 
glass phase should behave similarly to that in an ordered 
antiferromagnet (as the 9 held is frozen the spin auto¬ 
correlation is determined entirely by N), and will not 
show very striking glassy effects. In contrast there will 
be a wide range of relaxation times in a conventional spin 
glass. 

For PDW order the main proposed candidate to date 
are the cuprates, notably La2- x Ba x Cu20i. Close to 
x = g, there is an interesting window of tempera¬ 
tures between 4/\ and 16 K where the in-plane resistivity 
is immeasurably small while there is a non-zero c-axis 
resistivity^. Further the Meissner effect itself onsets only 
below 4 K. It has been suggested that this behavior may 
be explained by a PDW order pinned by impurities, and 
with frustrated c-axis Josephson coupling. Our results 
show that there are two possible fates of the PDW at 
weak disorder. In the PDW glass phase there are no 
frozen disorder-induced superconducting vortices. Con¬ 
sequently we expect a Meissner effect in this phase, and 
a non-zero critical current. In the fully disordered glass 
phase there are frozen random sign ^ vortices. This 
phase is then likely to behave similarly to a vortex glass. 
It will presumably have vanishing linear resistivity but 
a zero critical current and no Meissner effect. Thus it 
has the potential to explain the experiments within the 
PDW framework. An immediate consequence is the lo¬ 
cal breaking of time reversal symmetry at zero field due 
to the frozen vortices. It will be most interesting 
to look for this through scanning SQUID microscopy or 
other local probes of magnetism. 

Modulated superconductivity in the FFLO states, has 
of course, been discussed theoretically for decades. In 
recent years there have been suggestions of experimen¬ 
tal sightings of this state in two diff erent systems—first 
in the heavy fermion superconducto^EH CeCoIn$ and 
very recently in an organic superconductor®^. Both of 
these are very clean systems and hence our results on 
the effects of weak disorder may be directly applicable. 
CeCoIns is a 3D superconductor and therefore may be 
in a superconducting Bragg glass phase. The organic 
is quasi-two dimensional and thus will at best be in a 
phase with long range charge-4e superconducting order. 
In both materials it will be interesting to look for ^ flux 
quantization. 
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Appendix A: Functional Renormalization Group for 
the 3D Larkin-Ovchinnikov superconducting glass 


The starting point for the FRG analysis is the long- 
wavelength Hamiltonian 


Hfflo = ifcDW + H sc + H b (Al) 

Ecdw = K c /E^^+^E^-m ( A2 ) 


#sc = K s / Ew ™) 2 

Jx n 

H b = 2A / E • V<? 


(A3) 

(A4) 


where V{6) = V{—6) = V(9+n) is a symmetric, periodic 
function and /3 is the inverse temperature. In the SDW 
glass A = 0 and <f>, 9 decouple. The renormalization group 
equations for Hcdw were derived in Ref. [2D] as 


j f vw = ,vm + 


2V"(9) 2 -4V"(9)V"(0) 
&¥*! ’ 


(A5) 

(A6) 


where e = 4 — d and the tree level scaling of the elastic 
term has been absorbed into /?. In the present case of 
A ^ 0, it is clear that K s and A cannot be renormalized 
by V. The only modification to the RG equations is 
K c —>• K c = K c — X 2 /K s which can be absorbed by a 
simple rescaling so long as A 2 < K C K S . Thus the 3D 
FFLO glass exhibits the same universal properties as the 
3D SDW glass discussed in Sect. [IVJ 


Appendix B: Alternative lattice model 

In the main text we introduced the Hamiltonian 

H = — J Ni ■ Nj cos (9i — 9j + r/ij) — v cos(2 9i — ai ) 

<ij> i 

to capture the low-energy properties of the order parame¬ 
ters, in particular the structure of topological defects and 
the coupling to the disorder potential. A key feature of 
this Hamiltonian is a large redundacy corresponding to 
local gauge invariance under N r —► —Ni, 9i —> 9 t + 1 r. In 
some cases it is more convenient to adopt an alternative, 
equivalent formulation which makes this more explicit, 
and at the same time highlights the special role played 
by single dislocations. To this end we introduce a model 
in terms of N and bi = e l6i coupled to a Z 2 gauge field 
Oij with the Hamiltonian 


ffalt 


^ ) rr jj (^JsNi Nj -(- Jc cos {0j 9j 4- r/ij )^ 

ij 

—v E cos(2 9i — (Bl) 
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This has the same gauge invariance as H provided that 
Ni —> —Ni, Oi —> 8i + 7r is accompanied by u^- —► —<r,j 
for all sites j connected to i (the familiar ‘star’ trans¬ 
formation in Z 2 gauge theory). In the ground state one 
may choose the gauge cr,;,- = 1 to see that H and H a \ t 
yield the same energy for smooth fluctuations of N and 
e l6 . Moreover, H a \ t clearly allows 2tt vortices in N and 
double dislocations. 

In addition, there may be ‘visons’ i.e. ir flux configu¬ 
rations in the ct,; ; which cannot be removed by a gauge 
transformation. This flux is seen by both N and e lS , and 
therefore induces a single dislocation tied to a half-vortex 
in N. Thus the structure of topological defects in H and 
H a ^ is identical. Since the divergent contribution to the 
energy of any allowed defect is determined by the elastic 
terms in the Hamiltonian, H and H a \ t describe the same 
physics at long length scales. 

Appendix C: Estimates of length scales in 2D 

The presence or absence of isolated dislocations is de¬ 
termined by the balance between elastic energy cost and 
energy gain due to disorder. The elastic energy cost of a 
defect where 9 winds by 2Tvnn c and <j> winds by 2nm s 

Em c ,m s = n(K c ml + K s m 2 s ) log L/a. (Cl) 


To estimate the scale where defects proliferate, we drop 
the subleading factor and set p{£,m c ,m s ) = 1, obtaining 

c c J^l(^ 2 o+Ks/K c m^logi L /a 

U,m, =UeV 4Cm c (C7) 


Appendix D: Long distance physics for large -jf- in 2 d 

In this Appendix we sharpen the arguments justifying 
the absence of single dislocations at long wavelengths for 
large ^ in 2D. This is most conveniently done within 
the alternative formulation of the model introduced in 
App. 0 Although fJ a it can be fruitfully used to discuss 
all length scales/parameter regimes, we will use it here 
as an effective model for the large regime at length 

scales longer than ■* (the scale of doubled dislocation 
proliferation) discussed in the main text. Thus we will 
take the lattice spacing to be of order . The long 
length scale physics is then captured by this model with 
v —» 00 . Note that the core energy of a single dislocation 

must be taken to be E cole ~ K c In 
When t)-)oowe have 



The energy gain Vdis(®) due to disorder depends on the 
position x of the vortex and the particular disorder real¬ 
ization. The distribution of Vd; s for a constant variance 
of the random force is given bjEU 


Pm c (kdis) — f— exp 

CTTOcV ZTT 

(7 2 = 27 tD f K 2 log — 


vL A 

2a 2 m1) 


(C2) 

(C3) 


In the presence of random-fields, the variance of the ran¬ 
dom force is itself scale-dependent Dp(L) ss Clog^- 
where C is temperature dependent and > a is the 
Larkin length. In this case, a reasonable approximation 
consists of replacing the variance Dp of the random force 
by its average on a logarithmic scale 


Dp -A Dp(L) 


I 



Df{R) 

R 


dR 



L 


(C4) 


The probability p(L) that it a vortex exists in a volume 
Lr is given by 


P(L) 



Em c ,m s 


exp 



El 


2mlc 


(C5) 

(C6) 


CX-i 

T 



(Dl) 


with Si = d=l. Then the Hamiltonian becomes 


-^eff 


Gi i 

ij 


J s Ni ■ Nn + J c SiSj cos 




(D2) 

We define J Ci y = J c cos 2 ° J + . With rjij,ai ran¬ 

dom uncorrelated variables, J c ^j will have a probability 
distribution symmetrically distributed zero, and will be 
uncorrelated between different sites. In the absence of 
the <r l3 the second term describes a 2D Ising spin glass. 
In the configuration of J Ci ij there will with some probabil¬ 
ity be several frustrated plaquettes (where the product 
around the plaquette of J c ^j will be negative). If the 
gauge field a^ adjusts itself to “unfrustrate” those frus¬ 
trated plaquettes the Ising subsystem will gain energy 
AAising = A J c for some constant A. However such a 7r- 
flux of (Jij nucleates a 7r disinclination in N which costs 
( f (2) \ 

energy E corft ~ A c ln I I in addition to the elastic 


energy of distorting the N upto the typical separation 
between two such frustrated plaquettes. Thus the cost 
of nucleating a 7r-disclination in N overwhelms any en¬ 
ergy gain from unfrustrating the Ising spin Sj. Single 
dislocations are therefore suppressed in this regime. 
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